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ABSTRACT 

We determine the global star formation rate density at 0.7 < z < 1.9 using emission-line 
selected galaxies identified in Hubble Space Telescope Near Infrared Camera and Multi-Object 
Spectrograph (HST-NICMOS) grism spectroscopy observations. Observing in pure parallel 
mode throughout HST Cycles 12 and 13, our survey covers ~ 104 arcmin 2 from which we 
select 80 galaxies with likely redshifted Ha emission lines. In several cases, a somewhat weaker 
[OIII] doublet emission is also detected. The Ha luminosity range of the emission-line galaxy 
sample is 4.4 x 10 41 < L(Ha) < 1.5 x 10 43 erg s . In this range, the luminosity function is well 
described by a Schechter function with <p* = (4.24 ± 3.55) x 10" 3 Mpc -3 , L* = (2.88 ± 1.58) x 
10 42 erg s" 1 , and a = -1.39 ± 0.43. We derive a volume-averaged star formation rate density 
of 0.138 ± 0.058M Q yr _1 Mpc~ 3 at z = 1.4 without an extinction correction. Subdividing the 
redshift range, we find star formation rate densities of O.O88±O.O56M0yr _1 Mpc~ 3 atz= 1.1 and 
0.265 ±O.174M0yr _I Mpc~ 3 atz= 1.6. The overall star formation rate density is consistent with 
previous studies using Ha when the same average extinction correction is applied, confirming 
that the cosmic peak of star formation occurs at z > 1.5. 

Subject headings: galaxies: evolution - galaxies: luminosity function, mass function - galaxies: 
starburst - galaxies: high-redshift 



1. INTRODUCTION 

The evolution of cosmic star formation rate density (SFRD) and the stellar mass function are two key 
components required to describe galaxy evolution, representing current and past star formation activities 
respectively. Measurements of the SFR density at different redshifts (e.g., Madau et al. 1998; Steidel et al. 
1999; Arnouts et al. 2005; Schiminovich et al. 2005; Bouwens et al. 2006, 2007; Ly et al. 2007) suggest 
that the global SFR density peaks at 1 < z < 3. The study of the build-up of stellar mass density (e.g., 
Dickinson et al. 2003; Rudnick et al. 2003; Glazebrook et al. 2004; Fontana et al. 2004) also indicates that 

'Department of Physics & Astronomy, FPRD, Seoul National University, Seoul, Korea ; hjshim@astro.snu.ac.kr 
2 Spitzer Science Center, California Institute of Technology, Pasadena, CA 

3 Department of Physics & Astronomy, University of California, Los Angeles, CA 90095-1547; malkan@astro.ucla.edu 
4 Observatories of the Carnegie Institution of Washington, 813 Santa Barbara Street, Pasadena, CA 



-2- 



the redshift range 1 < z < 3 is the phase of massive galaxy formation, being the epoch of the strongest star 
formation. Most studies are in overall agreement that star formation decreases by a factor of 10 to 20 from 
z ~ 1 to z = 0. 

While the redshift range 1 < z < 2 is expected as the epoch of the strongest star formation, measure- 
ments of the SFR density over this redshift range are uncertain. First, most of the significant spectral features 
useful for redshift identification are in the near-infrared (0.7-2/mi) at z > 1, so few galaxies at z ~ 1-2 have 
spectroscopic redshifts. Second, the commonly used rest-frame ultraviolet or mid-infrared selections can 
be severely biased toward relatively unobscured or obscured star-forming galaxy populations. Therefore we 
performed a NIR spectroscopic survey using redshifted Ha emission lines at z > 1 , to select star-forming 
galaxies at the corresponding redshifts. 

Ha is known to be a robust measure of star formation, which is less affected by dust extinction com- 
pared to UV continuum (e.g., see review of Kennicutt 1998). At redshift 1 < z < 2, grism spectroscopy using 
Near Infrared Camera and Multi-Object Spectrograph (NICMOS) onboard Hubble Space Telescope (HST) 
offers a unique tool to sample Ha-selected star-forming galaxies. The first results of a NICMOS grism paral- 
lel survey were published using HST Cycle 7 data (McCarthy et al. 1999). Operating in pure parallel mode, 
the survey identified 33 emission-line galaxies over ~ 85 arcmin 2 of randomly selected fields. The Ha lu- 
minosity function at z = 0.7-1.9 was derived using the identified emission-line galaxies (Yan et al. 1999). 
Hopkins et al.(2000) investigate the faint-end of the Ha luminosity function in more detail by performing 
deeper pointed grism observations. This NICMOS parallel grism survey resumed after the installation of 
the NICMOS cryocooler and continued until 2005, including hundreds of observations throughout Cycles 
12 and 13. This extensive amount of new data enables the construction of a new, more robust Ha luminosity 
function at z = 1-2. 

In this paper, we present the Ha luminosity function at 0.7 < z < 1.9. Possible evolution of the lu- 
minosity function between redshift 1 and 2 is also investigated, by constructing luminosity functions at 
0.7 < z < 1.4 and 1.4 < z < 1.9 separately. The global star formation rate density inferred by the Ha lu- 
minosity function is compared with the values from previous studies, in context of the galaxy evolution. 
Throughout this paper, we use a cosmology with £Im = 0.3, = 0.7, and Hq = 71kms _1 Mpc _1 . 

2. DATA 
2.1. Observations 

All of the data presented in this study have been obtained using camera 3 of NICMOS onboard HST, 
taken in pure parallel mode with the G141 grism and the broad-band F110W/F160W filters. The original 
NICMOS camera 3 image is 256 x 256 pixel array with a pixel scale of 0.2" pixel -1 , providing a field of 
view of 51.2" x 51.2" (~ 0.75 arcmin 2 ). The observations were performed between October 2003 and July 
2004 (Cyclel2), July 2004 and January 2005 (Cyclel3). The fields are randomly selected, approximately 10' 
apart from the coordinates of the prime observation. The total exposure times for different fields vary from 
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768 seconds to 48000 seconds, while typical integration times ranged from ~ 2000 to ~ 30000 seconds. 
Every observed field has at least three dithered frames. To identify objects that are responsible for the 
spectra in slitless grism images, we obtained F160W (//-band) direct images before or after taking the grism 
images. For all Cycle 13 and several Cycle 12 fields, we also obtained F110W (/-band) direct images (For 
the information about the reduction and analysis of Fl 10W and F160W direct images, see Henry et al. 2007, 
2008). 

The G141 grism covers a wavelength range of 1.1 to 1.9/im, with mean dispersion of 8 x 10~ 3 //m 
pixel" 1 . The resolving power R is a function of the observing condition, including the variation of Point 
Spread Function (PSF) due to the changes in the optical telescope assembly and longer term changes in 
the internal structure (See McCarthy et al. 1999 for details). According to the previous observations, the 
nominal resolution is low: R ~ 100-200 (Noll et al. 2004). Thus, most of the lines are unresolved in this 
study. Assuming Ha emission line redshifted to 1 < z < 2 and R = 100 at A = 1 .5/xm, the smallest detectable 
rest-frame intrinsic equivalent width is 50-75 A, although there is some additional uncertainty resulting from 
the unknown relative strength of the [Nil] to the Ha line. 



2.2. Image Reduction 

The data reduction follows similar steps to previous studies using NICMOS grism data (e.g., McCarthy 
et al. 1999; Hopkins et al. 2000), including high background subtraction, one-dimensional spectra ex- 
traction, and wavelength/flux calibration. We start from the calibrated output *cal.fits files from the 
HST archive, which are corrected for bias and dark current removal, linearization, and cosmic ray rejection. 
Flat-fielding is not included in this stage but is done after the extraction of one-dimensional spectrum, since 
the flat-field strongly depends on the wavelengths in case of NICMOS grism. For images taken during the 
South Atlantic Anomaly, we apply SAA correction using saa_clear|_. After the SAA correction and 
the correction for any remaining differences in bias levels between each pedestal quarters, the main part 
of the data reduction is the removal of the high sky background at near-infrared wavelengths. We use two 
different methods to make the sky frame that will be subtracted, and select the better method for each case: 
i) median-combining all image frames taken and ii) median-combining only those image frames taken at 
close dates with the image frame that needs sky subtraction. For several fields taken during the period in 
which the sky background changes rapidly, the highly uneven background is not removed completely with 
the first method. In such cases, we use the second method to make a sky frame. At lease 9 frames taken at 
close dates were used to construct the sky frame, to prevent the increase in background uncertainty. Once 
the sky frame for each image is determined, the sky frame is subtracted from the observed images. Note that 
the construction of sky frames using image frames taken at close dates is newly introduced while previous 
studies used only one sky frame constructed by median-combination of all image frames (e.g., McCarthy et 
al. 1999). 



'http://www.stsci.edu/hst/nicmos/tools/post_SAA_tools.html 
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After the sky subtraction, we grouped the frames according to their 139 unique parallel fields and 
registered each group onto one plane. In order to measure the offsets between the dithered frames, we shift 
each frame by a series of Ax, Ay (±0. 1 pixel) shift values, subtract the shifted one from the reference frame, 
and find the shifts that minimize the difference through the iteration. The final shifts in x/y directions are 
less than 3 pixels in general. During this process, the bad pixels and hot pixels are masked out. For bad 
pixel masks, we combine the permanent bad pixel mask for the NICMOS camera 3 and the data quality 
flag image associated with the NICMOS raw data cube. Some "warm" pixels, which are missed in bad 
pixel masks are identified by eye and added in the final mask for correction. Final cleaned, sky-subtracted, 
coordinate-registered frames are added to construct the reduced mosaic two-dimensional image for each 
pointing. 

From the reduced two-dimensional spectra image, we extract one-dimensional spectrum and perform 
flux/wavelength calibration using the NICMOSlook softward^l developed by STEC-F (Freudling 1999). We 
first identify the position of the emission-line galaxy in the direct image (F160W) through visual inspection 
on the two-dimensional spectra image. From the interactively determined position and size of the object, we 
define the extraction aperture and the background region considering the offset between the grism and the 
direct image, and extract the spectrum. 

We apply a correction for the wavelength dependent pixel response to the extracted spectrum using an 
inverse sensitivity curve. For wavelength calibration, we first extract the spectrum of a bright point source 
in each grism image before extracting the spectrum of any emission-line galaxies. The bright point sources 
reproduce the significant cutoff in short/long limits (l/xm, 1.9/im) in inverse sensitivity curve with a high 
S/N, because stellar spectra are flat in the G141 bandpass. Final wavelength calibration for the spectrum 
of emission-line galaxy is done by adjusting slight offsets between the spectrum and the overlaid sensitivity 
curve. Note that sometimes the wavelength calibration varies for different locations in a camera field of 
view. The uncertainty in wavelength calibration due to the remaining distortion effect is about ~ 0.02/xm, 
i.e., in general there is a systematic redshift determination error of Az ~ 0.03. Also, there are a few cases 
where the emission-line object lies near the edge of the image. In these cases, the counterpart of the object 
is not found in the associated direct image, which causes a large uncertainty in wavelength calibration up 
to ~ O.l^m. The redshift uncertainty for an object lying at the image edge is Az ~ 0.15. Finally, the 
extracted one-dimensional spectrum is flux-calibrated by dividing the pixels values ([DN/s]) by the G141 
grism inverse sensitivity curve ([DN/s]/Jy). 

Our reduction method is comparable with that used in Hubble Legacy Archive (Freudling et al. 2008). 
The consistency between our spectra and the reduced spectra in Hubble Legacy Archivd_| confirms the 
existence of emission lines for the sample galaxies, although we find matches for only a few objects which 
are the most luminous. Most of our emission line galaxies are relatively faint, and are therefore missed by 
the less rigorous reduction in the Hubble Legacy Archive. 



2 http://www.stecf.org/instruments/NICMOSgrism/nicmoslook/nicmoslook/index.html 
3 http://hla.stecf.org 
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2.3. Area Coverage and Depth 

This work is an extension of the previous NICMOS grism survey for emission-line galaxies from the 
Cycle 7 parallel observation program (McCarthy et al. 1999), performed during Cycle 12 and Cycle 13. In 
Cycle 12, our survey targeted 130 different coordinates with different exposure times. From the original 130 
fields, 26 are excluded in the final analysis since the fields are either too crowded (i.e., the stellar densities are 
over ~ 50 arcmin -2 , M31/SMC fields), have high galactic foreground extinction (Taurus Molecular Cloud 
fields), or are damaged by latents produced by very bright objects observed just before the image exposure. 
In Cycle 13, we excluded 10 fields out of the 45 fields initially targeted for the same reasons. Therefore the 
total survey area is ~ 104 arcmin 2 , covering 139 different fields randomly distributed over the sky. 

We compare the exposure times of these 139 useable parallel fields taken during Cycles 12 and 13 with 
the earlier Cycle 7 survey in Figure [T^. Although the number of the deepest exposures is not significantly 
increased in Cycles 12 and 13, the number of total pointings is nearly double those from the Cycle 7 (Cycle 
7 data comprises 85 pointings over ~ 65 arcmin 2 ; McCarthy et al. 1999). In particular, the number of 
pointings with medium exposure times (2000-10000 seconds) have increased substantially. 

We illustrate the distribution of 5a line flux limits in Figure [Q). The rms noise in spectra-free regions 
is measured over a 4-pixel aperture for extracting one-dimensional spectrum. Thus, this "line flux limit" 
reflects the flux limit of a line added to the underlying continuum. As is expected from the exposure time 
comparison between Cycle 7 and Cycle 12/13 (Figure [1^), the line flux limits distribution of our data shows 
similar trend with that of Cycle 7. Though it is true that the fields with longer exposure time have fainter 
line flux limits, the depth of an image is not necessarily a simple function of the exposure time of an image. 
Instead, the flux limit is much more dependent on the flatness of the background, i.e., non-Poisson noise 
caused by imperfect background subtraction. Therefore, below 5a flux limits, we don't reliably identify 
emission lines because significant residuals remain from dark subtraction and other large data artifacts. 

3. EMISSION-LINE GALAXIES 
3.1. Identification 

We identify the emission-line galaxy candidates on the two-dimensional spectra image through visual 
inspection, prior to the extraction of a one-dimensional spectrum using NICMOS look (Section 2.2). We 
carefully compare the grism images and the direct images, identifying the galaxies in the direct images 
that are responsible for the emission-line in the spectra. We cross-check this method by identifying the 
same emission-line galaxies among multiple different authors. Some major obstacles in the identification 
of the emission-lines are the existence of zero-order images, the remaining background patterns, and the 
occasional image artifact. The zero-order image appears at a location ~ 27" apart from the end of the first 
order spectrum, so it can be recognized in most cases through the inspection of the first-order spectrum or 
the inspection of direct images. However in the middle portion of the detector, it is difficult to tell whether 
the point-like feature is a zero-order image or a strong emission-line with faint continuum. Figure [2] shows 
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two typical pairs of direct and grism images in our data and the identified emission-line galaxies. Zero-order 
images, the first and the second order spectrum, and the identified emission-line features are marked. 

Over the ~ 104 arcmin 2 surveyed in Cycle 12 and Cycle 13, we identify 80 emission-line galaxies. As 
we mentioned in Section 2.3, we only classify an emission-line as real if the line is significant at > 5c levels. 
TableQ]presents the coordinates, redshifts, line fluxes, observed equivalent widths, and available photometry 
of all 80 emission-line galaxies. The emission-line galaxies are distributed over 53 different fields, with 23 
fields containing more than one emission-line galaxy. 

The most likely candidates for the identified emission lines in our grism survey are hydrogen lines (Ha, 
H/3) and oxygen lines ([Oil] 3727A, [OIII] 5007A) redshifted to z > 0.7. Because of the broad wavelength 
coverage of G141 (1.1-1.9/um), we would expect to see both [OIII] and Ha in galaxies at 1.2 < z < 1.9 and 
both [Oil] and [OIII] in galaxies at 1.95 < z < 2.8. Only a single line is predicted for Ha at z < 1.2 and 
[Oil] at z > 2.8. Except for ten possible cases discussed in Section 3.3, we do not convincingly detect more 
than one emission line in most of the galaxies. That is, the derivation of redshift depends on only a single 
strong emission line in most cases. We believe that most of these single lines are Ha considering their large 
equivalent widths. 

We can estimate the possibilities of the emission lines being emission lines other than Ha using the 
expected equivalent widths of the lines, since our identification is limited to emission lines with EW (rest- 
frame) > 40-50A. First, the possibility of H/3 line is removed since the average equivalent widths for H/5 in 
star-forming galaxies are known to be relatively small (5-10A, Brinchmann et al. 2004). The next strongest 
line after Ha in terms of equivalent width is [Oil] 3727 A line. Identifying the emission line as [Oil] 3727 A 
requires the redshift of the object to be 1.95 < z, < 4. Considering the typical magnitude of our emission-line 
galaxies, the My should be ~ -24 mag if the galaxy is at z ~ 3. The estimated number of z ~ 3 galaxies 
with My < -24 over our survey volume is less than two according to the V-band luminosity function of 
z ~ 3 galaxies (Shapley et al. 2001), thus the possible contamination rate by [Oil] lines at z ~ 3 is less 
than 3 percents (2/80). Furthermore, most [Oil] contaminants will be in the 1.95 < z < 2.8 range, where 
[OIII] 5007 A line should also appear in our spectra. We suggest this possibility for one object, J0333 10.66- 
275221.4a presented in Figure@]in Section 3.3. This object is also included in TableQ]since the possibility of 
the line being [Oil] instead of Ha is still uncertain. Finally, the last possibility for contamination is the [OIII] 
5007 A line. Yet as we mentioned in the previous paragraph, [OIII] 5007 A emission line is accompanied 
with either Ha or [Oil] 3727 A at 1.2 < z < 2.8. Therefore the possibility of [OIII] being the single, strongest 
emission line of the galaxy is very low. 

Besides these strong candidates, lines with longer wavelengths are also possible, including Hell 10830A. 
However, these are not likely due to their weakness and the lack of accompanying nearby lines like [SIII] 
at 9069,9545A. The follow-up optical spectroscopy targeting 14 emission-line galaxies presented in Mc- 
Carthy et al.(1999) revealed that 9 of the emission-line galaxies are truly at z = 1-2, supporting the identified 
emission lines through NICMOS grism spectroscopy are mostly redshifted Ha (Hicks et al. 2002). The 
remaining 5 sources had no identifiable emission lines and were therefore unconfirmed, but not ruled out. In 
conclusion, we find that contamination from other emission lines is small (< 5%). Nonetheless, we include 
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this source of uncertainty in the derivation of the Ha luminosity function (Section 4). 

In addition to mis-identification of emission lines, there is also a possibility of contamination by Ha 
emission from AGN. We estimated the fraction of these contaminants using AGN luminosity function at 
similar redshifts and magnitudes. At z — 1, the number density of Type-1 AGN with (Mb) ~ -23 mag, 
corresponding to the median magnitude of our emission-line galaxies, is ~ 10~ 5 Mpc~ 3 (Bongiorno et al. 
2007). Our effective survey volume is ~ 3.14 x 10 5 Mpc 3 at this redshift range, which would indicate 
3^1 AGNs in our survey. This is actually the lower limit on the AGN number density in our data, since 
for a given //-band continuum magnitudes, AGN candidates are likely to be brighter and thus easier to be 
included in the grism survey sample compared to normal star-forming galaxies, due to the larger intrinsic 
equivalent width of the Ha line. Using the follow-up optical spectroscopy of the objects selected in the 
previous NICMOS grism survey, Hicks et al.(2002) suggested a new diagnostic for Seyfert 1 galaxies with 
L(Ha) and Ha equivalent width. We have five objects within the conservative cut of log(Ljj a )> 42.6 erg 
s -1 and EW(Ha)> 100A, while we have five more objects with EW(Ha)> 100 A and the luminosities of 
42.5 <log(L H J< 42.6 erg s" 1 . We assume these numbers represent the possible uncertainties caused by 
AGN contamination. 

To summarize, i) all the galaxies in the brightest bin (log(Ljj a )> 42.8 erg s -1 ) can be considered AGNs 
according to the criteria of Hicks et al.(2002). ii) ~ 40% of the galaxies in the bin of 42.5 <log(Ljj a )< 42.8 
erg s -1 are possible AGNs with large equivalent widths. These uncertainties are included in Table |2j and in 
the derivation of the Ha luminosity function. 

3.2. Ha-Derived Star Formation Rates 

Since Ha and [Nil] 65 83 A, 6548 A are not deblended at the resolution of NICMOS grism, we correct 
the derived Ha luminosity for [Nil] contribution when we use Ha luminosity as measures of star formation 
rate. The flux ratio used in correction is [Nil] 6583/Ha = 0.3 and [Nil] 6583/[NII] 6548 = 3 (Gallego et 
al. 1997). The Ha luminosities used throughout this paper are derived from this corrected Ha flux, derived 
using the following formula: / c (Ha) = 0.71 x /(Ha+ [Nil]). The star formation rates are derived from the 
corrected Ha luminosity using the formula of Kennicutt (1998), which assume a Salpeter IMF between 
0.1-100 Mq. Note that the star formation rate in TableQ]is not corrected for dust extinction. 

In Figure [3j we show the distribution of Ha-derived star formation rates as a function of redshift 
and Mr. The Ha-inferred star formation rates of the identified emission-line galaxies are 2-200 M yr _1 , 
comparable to or larger than the UV-estimated star formation rates of typical z ~ 3 Lyman break galaxies 
with relatively little extinction (Shapley et al. 2001). We do not apply any dust extinction correction for the 
star formation rate of individual galaxies or our derivation of Ha luminosity function, although (Ay) ~ 1 mag 
is expected for Ha-selected star-forming galaxies (e.g., Kennicutt 1992). We only use this Ha extinction 
correction for our estimates of star formation rate density evolution (see Section 5). 

The redshifts of our sample galaxies span the range of 0.7 < z, < 1.9, with a steep decrease at the low- 
redshift (z < 0.9) and high-redshift ends (z > 1.8). This is expected from the sharp end of the G141 grism 
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response curve. The median redshift is (z) = 1 .4. The inhomogeneous redshift distribution is considered in 
the derivation of the luminosity function, but the effect is small. 

Figure [3}) illustrates the median absolute magnitude of M R ~ -22.5 for our sample galaxies. The 
magnitude is comparable with Mg = -22.8 for galaxies at z = 1-1.2 (Ryan et al. 2007). Therefore our 
galaxies have typical luminosities around M*, across the entire redshift range. Despite considerable scatter, 
Mr can be used as stellar mass indicator for a galaxy. Thus, the higher SFR for the brighter galaxies may 
indicate that between 1 < z < 2, the star formation is higher in the more massive galaxies in our sample. 

3.3. Spectra of Possible Double-Line Objects 

We discussed in Section 3.1, the validity of identifying the single emission lines as Ha. If more than 
one emission line (i.e., two emission lines) exist, in nearly every case we identify them to be redshifted Ha 
and [OIII] 5007 A based on the wavelength ratio between two emission lines. Using the known [OIII]/Ha 
EW ratios of star-forming galaxies at z ~ 0.7, we predict the number of galaxies showing both Ha and 
[OIII]. For a typical [OIII]/Ha EW ratio of ~ 0.^0 and our detection limit of EW(emission line)> 40-50A, 
the Ha equivalent width required for [OIII] 5007 A line detection is > 120A. Additionally, in order to show 
both [OIII]/Ha lines in grism spectrum, the redshift of the object should be 1.2 < z < 1.9 which corresponds 
to ~ 60% of the total survey volume. Thus, the expected number of galaxies with both Ha and [OIII] 5007 A 
emission lines in this survey is ~ 11. 

The actual number of galaxies in our sample with possible [OIII] and Ha emission is eight or nine of 
the ten double-line objects. The line identification for one object is uncertain [J0333 10.66-27522 1.4a], and 
the other one [J121900.59+472830.9c] is likely a zero-order contaminant. The number is in good agreement 
with the expected ~ 1 1 Ha/[OIII] emitters. Figure |4] shows one-dimensional spectra of ten objects with 
plausible double-line features. 

J01410? '.80-652849. la - This object shows two emission lines at 1.73/um and 1.32/im although the 
line at 1.32^m is only marginally detected. These lines are likely redshifted Ha and [OIII] at z = 1.63. We 
could not derive a reliable flux for [OIII] line. According to the equivalent width and the luminosity of 
Ha line, this object is not considered as a Seyfert 1 galaxy (see Section 3. 1). The extended morphology 
(SExtractor CLASS_STAR < 0.8) also disfavors the possibility of this object being an AGN. 

J015240.44+005000.2a - This object shows two emission lines at 1.52/zm and 1.16/mi and both lines 
are significant. The lines are likely redshifted Ha and [OIII] at z = 1.31. The large Ha luminosity and the 
equivalent width, as well as the point-like morphology suggest that this object may be an AGN. The broad 
and asymmetric shape of the line at 1.16^m could be a blend of [OIII] 5007 A and H/3. 



4 This is the ratio calculated for the NICMOS grism emission-line objects using the follow-up optical spectroscopy (McCarthy 
et al. 1999; Hicks et al. 2002). The ratio is also comparable with the [OIII]/Ha EW ratio of emission-line galaxies in HST/ACS 
grism parallel survey (Drozdovsky et al. 2005) at 0.5 < z < 0.7. 
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J0333 10.66-27522 1 .4a - This object shows two emission lines at 1.32/um and 1.76/xm, but the iden- 
tification of these two lines is difficult because the wavelength ratio is uncertain. The measured flux ratio 
between the two lines EW\je im IEW\32fj,m is l ess than 1> so Ha would be weaker than [OIII]. On the other 
hand, the two lines may be explained more easily by [Oil] and [OIII] lines redshifted to z = 2.52. The mor- 
phology of this object is clearly extended, and the relatively smaller object size (radius of ~ 0.5") compared 
to other emission-line galaxies also suggests the possibility of this object at z > 2. The optical photometry of 
this object is available in MUSYC data (Gawiser et al. 2006), although the quality of SED fitting is low and 
the derivation of photometric redshift is difficult for this object. In Table [U we identify the line at 1.76/im 
as the Ha redshifted to z = 1-68. However, this object is an example of possible mis-identification in con- 
struction of emission-line galaxy sample, thus we include it in the uncertainties of the luminosity function 
in Section 4. 

J033310.66-275221.4b - This object shows two emission lines at 1.83/um and 1.39/um although the 
line at 1.39 fim is only marginally detected. The lines are likely redshifted Ha and [OIII] at z = 1.788. 
Judging from the clearly extended morphology and the combination of EW and luminosity of the Ha line, 
we conclude that this object is a star-forming galaxy rather than AGN. 

J121900.59+472830.9c - This object has a very bright broad line at 1.77/zm, and a weak line at 
1.31/mi. These may be redshifted Ha and H/3 at z = 1.7. However the line at 1.77/mi is strong (corre- 
sponding to a luminosity of L(Ha) = 1.8 x 10 43 ergs _1 , roughly ~ 5L*) and broad, so it might be from a 
zero-order image of an adjacent galaxy. The location of the potential counterpart for this possible zero-order 
image is out of the field of view of our direct image, so we cannot confirm clearly whether this is a zero-order 
image or a true emission line. On the other hand, if this line is indeed Ha, we measure a B aimer decrement 
of Ha/H/3 = 9.1, implying A v = 2.09. 

J122512.77+333425.1a - We see two lines at 1.57 fim and 1.19/um, which we identify as redshifted Ha 
and [OIII] at z = 1.39. The measured flux ratio between the two lines is f\.$ijj, m lf\.\9iim ~ 1-06, which falls 
within the possible range of Ha/[OIII] values for star-forming galaxies. Moreover the equivalent width and 
the luminosity of Ha line remove the possibility of this object being a Seyfert 1 galaxy. The object shows 
relatively symmetric but not compact morphology in the direct (F160W) image, thus providing support that 
this object is a star-forming galaxy rather than AGN. 

J123356.44+091758.4a - This object shows two emission lines at 1.84/zm and 1.38/mi although the 
line at 1.38/im is very weakly detected. The lines are likely redshifted Ha and [OIII] at z = 1.8. The large 
size of this object, large Ha EW and large Ha luminosity suggest this object is a large galaxy being powered 
by AGN. 

J 125424.89+270147. 2a - This object shows two emission lines at 1.88/im and 1.43^m. The lines 
are likely redshifted Ha and [OIII] at z = 1.861. Judging from the clearly extended morphology and the 
combination of EW and luminosity of Ha lines, we conclude that this object is a star-forming galaxy rather 
than AGN. 

J161349.94+655049.1a - In addition to the significant line at 1.9/zm, the spectrum of this object shows 
a weak line at 1.43 fim. These lines are likely redshifted Ha and [OIII] at z = 1.9, although the line at 
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1.43/xm is only marginally detected. The line at 1.9/um has relatively large equivalent width (473A). The 
spatial scale of the object is relatively large and shows a sign of substructure, so we classify it as a bright 
star-forming galaxy. 

J213717.00+125218.5a - This object shows two lines at 1.68/um and 1.28/um, which are probably 
redshifted Ha and [OIII] at z = 1.56. The flux ratio between the two lines is /1.68/W 'fi.2S/jm ~ 1-55, although 
the [OIII] line is weak. The morphology of this object in the F160W image is not classified as a point 
source (with CLASS_STAR less than 0.9 in SExtractor output), although it does not show any clear signs of 
interaction. 



We use the 1/V max method (Schmidt 1968) in order to construct the Ha luminosity function. Although 
the redshift distribution of the galaxies is not homogeneous due to the sharp cutoff at low/high redshift end, 
we find the effect is negligible in the calculation of V max . We begin by calculating the maximum comoving 
volume V max over which each galaxy could lie and be detected, including corrections for all of our sample 
selection biases: redshift, flux, and the location in the image. The equation used is as follows (see eqn [1] 
of Yanetal. 1999) : 



In the above equation, the differential comoving volume dV / dz is integrated over the redshift range of 
[zl, z2]. The integration range is defined as zl = max(z/,z mm ), z2 = min(z/, , Zmax) while zi and Zh indicate 
the lower/higher redshift ends limited by the spectral cutoff of the G141 response curve (zi = 0.67 and 
Zh = 1-9). The quantities z m i n and Zmax are the minimum/maximum redshift the object can be detected. 
While max(z/,z m j n ) is zi in general, the maximum redshift the object can be detected is determined by the 
measured object flux, according to the equation Di(zmax) = Di(z)[fn a / fii m Y^ 2 where Di is the luminosity 
distance at redshift z, fn a is the measured flux of the identified Ha line, and f^ m is the limiting line flux 
of the image (see Section 2.3.). The spectral cutoff zi and zh vary according to the location of the object in 
NICMOS grism field of view. Covering the whole wavelength range of 1.1-1.9/im is only possible when 
the object lies at central portion of the image. Therefore, we used different zi and zh according to the object 
location in the grism field of view. 

The observed luminosity function needs to be corrected for sample selection biases, including redshift 
limits and the incompleteness. R(X) is the G141 inverse sensitivity curve, which corrects for the effect of 
the sharp spectral cutoff at the low- and high-redshift ends. Since we are using images with different depths, 
different incompleteness corrections must be applied for individual objects. C(/jj a ,/ij m ) is a factor for 
incompleteness correction assigned to each object, as a function of the observed Ha line flux and the flux 



4. Ha LUMINOSITY FUNCTION 




(1) 
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limit of the image. To estimate this correction factor, we first select several high S/N emission-line galaxies 
with different equivalent widths and wavelengths. After cutting out two-dimensional spectral templates of 
the selected emission-line galaxies, we dim the images by various factors to generate artificial emission-line 
galaxies with various line fluxes. The dimmed two-dimensional spectra are added to original NICMOS 
grism images at random locations, then we perform the one-dimensional spectral extraction, applying the 
same method that we described in Section 2.2. We repeat these steps for images with different /jj m to 
measure C(/jj Q ,/jj m ) as a function of galaxy location. We plot both uncorrected and corrected luminosity 
function in Figure [5^ in order to show the significance of the incompleteness correction. Finally, f2 is the 
solid angle covered in this survey. Q is also a function of /jj since our objects are gathered from multiple 
images with widely differing line detection depths. 

The source density in a specific luminosity bin of width A(logL) centered on the luminosity logL, is the 
sum of (1/Vmax) of all sources within the luminosity bin (i.e., logL,+i - logL,- = 2 x A(logL)). The variances 
are computed by summing the squares of the inverse volumes, thus the luminosity values and the error-bar 
in each bin is evaluated using the following equations. 



^(logL^— Yl T7~ — W 
A(logL) , , , ^max 
|logL-logz.,|<A(logz.) 



1 

0> : 



AlogL \ ^ 

\ |logL-logL,|<A(logL) 



max 



After calculating the size of the error-bars using the equation, we added additional uncertainties that 
might result from AGN contamination and line mis-identification (please refer to Section 3.1 for details). 
Note that the uncertainties in the luminosity function and the derived SFR density due to the large scale 
structure are less than ~ 2% (Trenti & Stiavelli 2008), given the large number of independent fields used in 
our survey. The derived luminosity functions from the whole sample (0.7 < z < 1 .9) is presented in Figured 
and Table|2j while the errors indicated already include additional uncertainties from AGN contamination and 
line mis-identification. The points from the previous NICMOS grism selected Ha emitters (Yan et al. 1999; 
Hopkins et al. 2000) are overplotted for comparison, after accounting for the cosmology differences. Over 
the luminosity range of 41.6 < logL(Ha)[erg s -1 ] < 43.5, our LF is consistent with those of the previous 
studies. The solid line in Figure [5^ is the best-fit Schechter LF to our data points, yet the line still fits 
points from other studies as well within the error bars. The plot shows a clear evolution of the Ha LF from 
z = 1.4 to the local Ha LF (dotted line), a result already well-known from the luminosity evolution in other 
wavelengths like the UV (e.g., Arnouts et al. 2005). 

The Ha luminosity function we have derived is for galaxies with EW (rest-frame) > 40A, as our survey 
is unable to confidently detect lines with lower equivalent widths. However, we note that our observed range 
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of EWs is comparable to that of z ~ 2 star-forming galaxies observed in Erb et al.(2006). This suggests that 
there may be substantial overlap between our sample and the one reported by Erb et al.(2006), although 
our sample likely contains some dust obscured galaxies that would be missed by the UV selection of Erb et 
al.(2006). 

The best-fit Schechter function parameters (Table [3]) are derived using the MPFIT packaged, which 
provides a robust non-linear least square curve fitting (e.g., Ly et al. 2007). The errors in Table [3] correspond 
to la uncertainty for each Schechter function parameter, and are derived using a Monte Carlo simulation. 
We generated a large number (~ 10000) of Monte Carlo realizations of our LF, with different [logL, log0] 
sets perturbed according to the uncertainties. Then we repeat the fit to find the best-fit parameters for each 
realization of the LE We find the faint-end slope is largely unconstrained by our data: a = -1.39 ±0.43. 
This is consistent with both the local Ha LF (a = -1.35; Gallego et al. 1995), and the deep NICMOS grism 
Ha survey (a = -1.86 ±0.14; Hopkins et al. 2000). In addition to the Schechter function fitting with varying 
a, we derived L* and (f>* with a being fixed to -1.39,-1.0 and -1.8 (Table© to cover all the possible range 
of a, and to investigate the effect of varying a on the total SFR density derived. 

Since the number of our sample galaxies (80) is more than twice of that from the previous studies, 
we can also test the evolution of Ha LF as a function of redshift between 0.7 and 1.9. We divided the 
sample galaxies into two redshift bins (0.7 < z < 1.4, 1.4 < z < 1.9), and derived the LFs separately. The 
LFs for two sub-samples at different redshift bin are shown in Figure [5t>. The luminosity function in the 
lower-redshift bin (0.7 < z < 1.4; (z) = 1.1) is plotted as triangles, while the luminosity function for the 
higher-redshift bin (1.4 < z < 1-9; (z) = 1.6) is plotted as squares. We see that there are more Ha-luminous 
galaxies in the higher-redshift bin than in the lower-redshift bin, which implies that L* , (jf , or both are larger 
at higher redshift. 

The LF values and the best-fit Schechter parameters for these two redshift bins are also listed in Table 
|2] and [3] We also show the derivation of the Schechter parameters for the faint-end slope fixed at a = 
-1.39,-1.0, and- 1.8 as this quantity is more difficult to constrain for these smaller sub-samples. The errors 
in cj)* and L* are derived by the same Monte Carlo method described above. For the cases where a is fixed, 
the uncertainties are artificially decreased, so we adopt the larger uncertainties derived when a is free. In 
the inset plot of Figure [5t>, we illustrated how <fi* and L* evolve from 1.4 < z < 1.9 to 0.7 < z < 1.4. The 
contours represent the la uncertainty range for the parameters (for a free and fixed, dot-dashed/solid line). 
The contours for parameters at 0.7 < z < 1.9 is also illustrated. 



5. EVOLUTION OF SFR DENSITY at 1< z < 2 

From the derived LF parameters, we evaluate the Ha-inferred star formation rate density at 0.7 < 
z < 1.4 (lower-redshift bin), and 1.4 < z < 1.9 (higher-redshift bin). The SFR densities are also listed in 
Table [3] The total integrated Ha luminosity density is calculated as L, ot = <j)*L*T(a + 2), then converted to 



3 http://cow.physics. wisc.edu/~craigm/idl/fitting.html 
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SFR density using the relation of Kennicutt (1998) : SF/?(//a)[M yr 1 ] = 7.9 x 10" 42 L(//Q)[ergs -1 ]. The 
relation assumes a Salpeter IMF between 0.1-100 M . The uncertainties in Table[3]are also derived through 
Monte Carlo method, by making a large realization of [<fi*, L*, a] sets, calculating L tot , and estimating 
la uncertainty from the distribution of L tot . As mentioned in previous section, the uncertainties on L* are 
underestimated for the cases where a is fixed. 

We compare our SFR density estimates with the results of other studies in Figure [6] The SFR points 
are drawn from individual references (Gallego et al. 1995; Yan et al. 1999; Hopkins et al. 2000; Moorwood 
et al. 2000; Perez-Gonzalez et al. 2003; Sullivan et al. 2000; Tresse & Maddox 1998; Tresse et al. 2002; 
Pascual et al. 2001; Fujita et al. 2003; Nakamura et al. 2004; Hippelein et al. 2003; Glazebrook et al. 1999, 
2004; Doherty et al. 2006; Ly et al. 2007; Shioya et al. 2008; Villar et al. 2008; Reddy et al. 2008; see 
compilation of Hopkins 2004), and were corrected to fit different cosmology. For extinction, we used the 
values given by the authors. The amount of the applied extinction collection is different for different studies, 
from A(Ha) ~ 0.3 mag to A(Ha) ~ 1.2 mag. This difference up to 1 magnitude produce the uncertainty in 
the SFR density points of a factor of ~ 2. If the authors do not provide information about the extinction, we 
apply the average extinction of (Ay) = 1 mag, i.e., (A(Ha)) ~ 0.85 mag (e.g., Kennicutt 1992), the factor 
widely adopted in previous studies (e.g., Hopkins 2004; Doherty et al. 2006). The extinction collection 
applied to our SFR density points (stars in Figure© is also (A(Ha)) ~ 0.85 mag. 

Our estimate of the volume-averaged SFR density at 0.7 < z < 1.9 is consistent with that of previous 
NICMOS grism studies over the same redshift range (points at z = 1.3 from Yan et al. 1999; z = 1.25 from 
Hopkins et al. 2000). Note that the integration ranges for Ha luminosity density are different from study to 
study - our study accepts L in the range [0, oo] while Hopkins et al.(2000) have integrated the LF over the 
range of 10 37 < L[erg] < 10 47 . Since the derived faint-end slope in our study is relatively flat (1.39 ±0.43), 
this difference of integration range makes little difference in the final SFR density value. If we restrict the 
integration range to 10 37 < L[erg] < 10 47 as in Hopkins et al.(2000), our result is decreased by only ~ 0.1%. 

Our points clearly place the peak epoch of the relatively "unobscured" star formation at z = 1-2. More- 
over, by examining the two points at z = 1.1 and z= 1.6, we can conclude that the peak of star formation 
must have occurred at redshifts higher than z = 1.5. This result is fairly robust, as it uses two sets of galaxies 
all identified with the same selection method. 



6. SUMMARY 

We have designed and executed a NICMOS grism survey, exploring the rest-frame optical universe at 
0.7 < z < 1.9. Through this program, we have identified a unique sample of emission-line galaxies over a 
significant cosmic volume, which enable us to study the relatively bright part of the Ha luminosity function 
at 0.7 <z< 1.9. 

Using Cycle 12 and Cycle 13 data, we probe ~ 104 arcmin 2 area, at 139 different locations throughout 
the sky. This corresponds to an effective comoving volume of ~ 3.14 x 10 5 Mpc 3 , almost two times larger 
than that of our previous observations (McCarthy et al. 1999). We identified 80 probable emission-line 
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galaxies, down to L(Ha) ~ 4.4 x 10 ergs" . Most of the emission lines are thought to be redshifted Ha, 
and their //-band magnitude distribution suggests that the identified emission-line galaxies are relatively 
bright (M ~ M*) star-forming galaxies at 0.7 < z < 1.9. 

We construct the Ha luminosity function from the line fluxes and the redshifts of these galaxies. From 
the integration of the luminosity function, the luminosity density and the star formation rate density are 
derived. We divide our sample into two redshift bins, one at 0.7 < z < 1.4 and the other at 1.4 < z < 1.9. 
The volume-averaged star formation rate densities at these two different redshift range are evaluated to be 
0.088 ± O.O56M yr _1 Mpc" 3 and 0.265 ± 0. 174M yr _1 Mpc" 3 respectively. The results are consistent with 
other Ha-derived SFR densities at similar redshifts. Using our unique sample, all selected by the same 
method, we find that the cosmic star formation history probed by these Ha measurements places the peak 
of star formation at z > 1.5, with a decrease in global SFR density from z=1.6toz=l.l. Although there 
remain uncertainties in the relative extinction, our study places firm constraints on the cosmic star formation 
history. 
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Table 1. Emission-line objects from the NICMOS grism survey 
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J002638.02+170117.9 


a 


00:26:38.32+17:01:20.7 


1.206 


2.653 


320 


23.13 


22.52 


42.24 


13.8 




b 


00:26:38.35+17:01:24.8 


1.300 


1.299 


230 


23.48 


22.61 


41.97 


7.4 


J013433.61+311818.2 


a 


01:34:34.38+31:18:05.3 


1.283 


3.372 
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42.38 


19.0 




b 
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a 
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_ 
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J013655.80+153720.4 


a 
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12.9 
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a 
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42.38 


18.8 
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b 
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a 
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a 
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b 
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a 
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J033309.36-275242.1 


a 
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_ 
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5.5 




b 
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1.779 
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_ 
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42.22 


13.0 
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a 


03:33:07.68 -27:51:47.1 
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_ 
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41.94 
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b 
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_ 
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42.20 
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a 
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1.465 
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325 
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21.45 


42.55 


28.2 
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a 


08:19:55.53+42:17:50.5 


1.122 
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20.82 


42.28 


15.0 




b 


10:47:02.36+12:29:35.6 


1.482 


3.175 


542 




21.79 


42.43 


21.3 


J104849. 15+462852.5 


a 


10:48:48.37 +46:28:57.0 


1.632 


3.292 


445 




23.49 


42.50 


24.9 


J111942.69+513839.1 


a 


11:19:40.61 +51:38:55.8 


1.733 


1.689 


219 


23.09 


22.00 


42.24 


13.8 


Jl 12413.26-170215.6 


a 


11:24:13.50 -17:02:14.1 


1.744 


0.719 


185 


23.94 


23.47 


41.87 


5.9 




b 


11:24:13.70 -17:02:33.2 


1.349 


2.517 


443 


23.53 


22.34 


42.28 


15.1 



Table 1 — Continued 



Field Object" a (J2000) <5 (J2000) redsbifW Flux/; a b W obs F110W c F160W c logL Ha d SFR d 

(1CT 16 erg s" 1 cm 4 ) (A) (mag) (mag) (erg s -1 ) (Mq yr~' ) 



J112414.50-170137.3 


a 


11 


:24 


:15.58 -17:01:45.9 


1.667 


1.041 


185 


23 


.79 


23.02 


42.01 


8.1 


Jl 12846.94+641459.7 


a 


11 


:28 


:44. 19 +64:15:06.3 


1.005 


1.467 


137 


23 


.62 


23.75 


41.89 


6.1 




b 


1 1 


:28 


:48.60 +64:15:00.5 


1.559 


2.655 


177 


22 


.58 


23.33 


42.38 


19.0 


J 1205 13. 7 1-073 124.0 


a 


12 


:05 


:13.33 -07:30:54.9 


1.550 


0.722 


190 






25.18 


41.81 


5.1 


J121900.59+472830.9 


a e 


12 


:18 


:58.73 +47:28:16.6 


1.679 


1.751 


236 








42.24 


13.7 




b 


12 


:18 


:59.24 +47:28:10.9 


1.462 


4.513 


1 15 


22 


.92 


22.16 


42.58 


29.8 




c 


12 


:18 


:59.48 +47:28:10.8 


0.996 


3.779 


200 


21 


.63 


21.06 


42.29 


15.5 


J121936.36+471950.6 


a 


12 


:19 


:35.78 +47:19:35.8 


1.186 


1.391 


220 






24.56 


41.95 


7.1 


J122226.64+042914.2 


a 


12 


:22 


:25.25 +04:28:35.9 


1.601 


0.556 


106 






22.88 


41.72 


4.1 


J122246.79+155704.0 


a 


12 


:22 


:47. 16 +15:56:45.2 


1.697 


0.853 


126 


22 


.94 


22.09 


41.93 


6.8 


J122512.77+333425.1 


a 


12 


:25 


:11. 12+33:34:31.7 


1.392 


1.950 


148 


22 


.68 


21.72 


42.19 


12.1 




b 


12 


:25 


:13.26 +33:34:14.3 


1.479 


1.402 


209 


23 


.89 


22.55 


42.08 


9.4 


J122908.41+015420.1 


a 


12 


:29 


:08. 56 +01:54:08.5 


1.752 


1.868 


171 






20.96 


42.29 


15.4 


J123041.86+121509.1 


a 


12 


:30 


:42.28 +12:15:26.9 


1.475 


2.437 


188 






22.40 


42.31 


16.3 


J123356.44+091758.4 


a 


12 


:33 


:56.24 +09:17:54.0 


1.800 


14.161 


371 






20.75 


43.18 


120.9 


J124339.45-341843.1 


a 


12 


:43 


:40.88 -34:18:20.6 


1.050 


3.115 


146 


22 


.58 


21.87 


42.24 


13.7 


J125424.89+270147.2 


a 


12 


;54 


:24.63 +27:02:05.2 


1.861 


2.837 


211 






22.51 


42.50 


25.2 


J132745.26-31 1537.8 


a 


13 


:27 


:44.50 -31:16:07.5 


1.240 


2.566 


251 


23 


.56 


23.13 


42.24 


13.8 


J132820.08-3 13744.8 


a 


13 


:28 


: 19.70 -31:37:42.2 


1.280 


1.474 


274 


22 


.68 


21.84 


42.02 


8.3 


J135835.94+623046.5 


a 


13 


:58 


:38.24 +62:30:55.8 


1.476 


0.638 


203 


23 


.44 


22.88 


41.73 


4.3 


J140214.22-1 13634.5 


a 


14 


:02 


:13.67 -11:36:53.1 


1.514 


0.711 


215 


23 


.09 


22.85 


41.79 


4.9 


J141833.44+250745.0 


a 


14 


:18 


:33.84 +25:07:45.1 


1.419 


1.629 


215 


24 


.48 


24.39 


42.12 


10.4 


J161349.94+655049.1 


a 


16 


:13 


:52. 19 +65:50:51.4 


1.904 


6.688 


473 






23.66 


42.89 


61.2 




b e 


16 


:13 


:47.61 +65:50:35.0 


1.748 


0.692 


195 








41.86 


5.7 


J175907.13+664454.3 


a 


17 


:59 


:03.06 +66:45:13.0 


1.331 


1.620 


178 


24 


.56 


24.72 


42.08 


9.5 




b 


17 


:59 


:09.47 +66:44:48.9 


1.614 


1.618 


190 


25 


.12 


24.13 


42.18 


12.1 




c 


17 


:59 


:1 1.37 +66:44:41.9 


1.446 


0.640 


131 


24 


.11 


23.67 


41.72 


4.2 


J213717.00+125218.5 


a 


21 


:37 


:18.38 +12:51:34.8 


1.561 


6.620 


226 






20.85 


42.78 


47.4 


J220239.22+185112.4 


a 


22 


:02 


:38.87 +18:51:37.9 


1.406 


2.759 


312 






23.40 


42.34 


17.4 




b 


22 


:02 


:39.09 +18:51:29.3 


1.505 


4.187 


484 


23 


.72 


23.13 


42.56 


28.7 




c 


22 


:02 


:40.03 +18:51:24.5 


1.663 


1.678 


242 


24 


.62 


23.60 


42.22 


13.0 


J222640.90-722919.8 


a 


22 


:26 


:41.14 -72:29:14.3 


1.692 


1.477 


191 






20.95 


42.17 


11.7 


J225916.39-345408.5 


a 


22 


:59 


:16.51 -34:54:11.5 


1.507 


2.931 


254 






22.39 


42.41 


20.1 


J230342.56+085617.2 


a 


23 


:03 


:42.20 +08:56:42.9 


1.355 


5.527 


156 


22 


.46 


21.59 


42.62 


33.3 




b 


23 


:03 


:43.68 +08:56:07.4 


1.164 


0.799 


581 


23 


.33 


23.00 


41.70 


4.0 



a We assign the suffix a, b, or c to identify different objects in one field. 

b Fluxtf Q is the emission line flux. Since the [Nil] and Ha lines are not resolved in the resolution of NICMOS grism, Flux« Q 
represents f(Ha + [Nil]). 

C 'F1 10W/F160W magnitudes are in AB magnitudes, MAG_AUTO from SExtractor (total magnitudes of the galaxies). The conver- 
sion between AB and Vega magnitudes are F110Wy egl7 = Fl 10W^ B -0.73; VlbWNvega = F160Wab-1.31 (derived using NICMOS 
zeropoints in Vega magnitude system at http://www.stsci.edu/hst/nicmos/performance/photometry/postncs_keywords.html I. 

d Ha luminosity and the derived star formation rate are corrected for possible [Nil] contamination in the measured Ha line flux: 
f(Ha) = 0.71 X f(Ha + [Nil]) 

c The objects that do not have either Fl 10W/F160W photometry are objects lying near the image edge. Due to the locations, we 
could not find these objects in the direct image despite of clear emission lines (not thought to be zero order) for these galaxies. All 
galaxies without Fl 10W photometry do not have the corresponding Fl 10W images, i.e., there is no Fl lOW-dropouts or galaxies with 



very red (F110W-F160W) colors. 
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Table 2. The derived Ha luminosity function in grism survey 





z ~ 0.7-1.9 






z ~ 0.7-1.4 






z ~ 1.4-1.9 




logL 


<f> (10" 3 Mpc- 3 ) 




logL 


4> (10" 3 Mpc- 3 ) 


Ngal 


logL 


(f) (10" 3 Mpc- 3 ) 


Ngal 


41.8 


5.7115±2.8134 


20 


41.9 


3.6746 ±1.3633 


13 


41.85 


7.4484 ±4.6620 


13 


42.1 


3.2962 ±0.7774 


26 


42.3 


1.6812 ±0.4203 


17 


42.25 


3.3088 ±0.8610 


20 


42.4 


2.0319 ±0.5149 


21 


42.7 


0.3485 ±0.2248 


4 


42.65 


0.9663 ±0.3960 


10 


42.7 


0.6172 ±0.2522 


10 


43.1 


0.0871 ±0.0871 


1 


43.05 


0.1868 ±0.1868 


2 


43.0 


0.1202 ±0.1202 


2 














43.3 


0.0704 ± 0.0704 


1 


















80 






35 






45 



Note. — The uncertainties include Poisson errors, possible AGN contamination, and line mis-identification. 




Fig. 1. — (Top): The distribution of exposure times of the observed fields in Cycles 7, 12, and 13. The 
shaded histogram with dotted lines is the data presented in McCarthy et al. (1999). The open histogram 
with dashed/solid line is the accumulated distribution when Cycles 12/13 data are added to the existing 
data. That is, the blue dashed line indicates the sum of Cycle 7 and Cycle 12 data, while the red solid line 
indicates the sum of Cycle 7, Cycle 12, and Cycle 13 data. (Bottom) : The histogram of 5a line flux limits, 
within a 4-pixel aperture of the NICMOS grism data. Each pointing covers 0.75 arcmin 2 . The line flux 
limits are compared to those from McCarthy et al. (1999), which is drawn with a dashed line. 
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Table 3. Ha luminosity function parameters for emission-line galaxies 



z range N gat 


0*(lO" 3 Mpc- 3 ) 


a 


logL* 


SFRD (Moyr^Mpc" 3 ) 









0.788 ±0.881 


-1.78 ±0.46 


42.72 ±0.23 


0.088 ±0.056 


0.7-1.4 


35 


2.704 ±0.841 


-1.39 (fixed) 


42.48 ±0.12 


0.094 ±0.023 






5.046 ±1.125 


-1.0 (fixed) 


42.27 ±0.09 


0.074 ±0.014 






0.848 ±0.289 


-1.8 (fixed) 


42.79 ±0.12 


0.187 ±0.033 






1.745 ±1.578 


-1.90 ±0.39 


42.72 ±0.20 


0.265 ±0.174 


1.4-1.9 


45 


4.096 ±1.129 


-1.39 (fixed) 


42.54 ±0.10 


0.164 ±0.034 






7.106 ± 1.851 


-1.0 (fixed) 


42.37 ±0.09 


0.129±0.019 






1.491 ±0.495 


-1.8 (fixed) 


42.80±0.12 


0.334 ±0.063 






4.241 ±3.553 


-1.39 ±0.43 


42.46 ±0.19 


0.138 ±0.058 


0.7-1.9 


80 


3.303 ±0.512 


-1.39 (fixed) 


42.54 ±0.06 


0.132±0.012 






5.888±0.514 


-1.0 (fixed) 


42.38 ±0.03 


0.1 12 ±0.009 






1.122±0.188 


-1.8 (fixed) 


42.82 ±0.04 


0.269 ±0.021 



Note. — The values are calculated for £l m = 0.3, = 0.7, and h = 71kms ^pc . The 
adopted form of the luminosity function is cf)(L) = 4>*(L/L*) 1+a e\p(-L/L*). The SFR den- 
sity column is inferred from the Ha luminosity density, using S7 7 /?(//a)(M yr~ 1 ) = 7.9 x 
10~ 42 L(//a)ergs _1 from Kennicutt (1998). Here, the total luminosity density is derived using 
the luminosity function parameters : L t0 , = 0*L*T(a+2). 
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Fig. 2. — Typical pairs of direct (F160W) and grism (G141) images of our NICMOS parallel survey. Each 
field is 51.2" x 51.2". In the grism two-dimensional image on the right, we mark the zero-order image (Z), 
the first and second order spectrum (1st order / 2nd order), and the emission lines. The objects producing 
the emission lines are also marked as A, B, or C in the direct image. 
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Fig. 3. — {Left): Star formation rates of the emission-line galaxies as a function of redshift. The identified 
emission-line galaxies are distributed over 0.7 < z < 1.9, with no significant redshift peak. (Right): Star 
formation rates of the emission-line galaxies as a function of absolute magnitudes Mr. Mr magnitudes are 
derived from the observed //-band (F160W) magnitudes. 
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Fig. 4. — One-dimensional spectra of the galaxies which appear to have more than one emission line. In 
most cases, the emission lines are redshifted Ha and [OIII] 5007 A (see Section 3.3 for details). 
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Fig. 5. — (Left): Ha luminosity function over 0.7 < z < 1 .9 derived from our study of emission-line galaxies 
in NICMOS grism survey. The observed values are plotted as open circles, and the values after the incom- 
pleteness correction are plotted as filled circles. Error bars indicate Poisson (solid) and additional(<io?te<i; 
uncertainties from line mis-identification, contamination) errors. Compared are Yan et al.(1999), Hopkins et 
al.(2000) points also from NICMOS grism studies (cross/triangle respectively). The solid line is the best-fit 
Schechter luminosity function to our derived luminosity function points, while the dotted line indicates the 
Ha luminosity function of local galaxies (Gallego et al. 1995). (Right): Ha luminosity function derived 
at two different redshift range (0.7 < z < 1.4, 1.4 < z < 1.9). Overplotted solid and dashed lines are the 
best-fit Schechter luminosity function with the faint-end slope a fixed to -1.39. In the inset plot, the best-fit 
Schechter LF parameters for the two redshift-bin sub-samples are shown in addition to the Schechter param- 
eters for 0.7 < z < 1.9 samples. The solid contours indicate la uncertainties in L* and cf>* when faint-end 
slope a is fixed, while dot-dashed contours indicate la uncertainties when a is a free parameter. 
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Fig. 6. — The evolution of SFR density as a function of redshift. Our points are plotted as two filled stars 
at z = 1.1 and 1.6, produced from the integrated luminosity functions of the two sub-samples covering these 
redshift ranges. Also, the SFR density evaluated from the luminosity function over the entire range of 
0.7 < z < 1.9 is plotted at z = 1.4 (open star). All other plotted points are based on the Ha-derived SFR 
density from spectroscopy or narrow-band imaging, (e.g., Gallego et al. 1995; Yan et al. 1999; Hopkins et 
al. 2000; Moorwood et al. 2000; Perez-Gonzalez et al. 2003; Sullivan et al. 2000; Tresse & Maddox 1998; 
Tresse et al. 2002; Pascual et al. 2001; Fujita et al. 2003; Nakamura et al. 2004; Hippelein et al. 2003; 
Glazebrook et al. 1999, 2004; Doherty et al. 2006; Ly et al. 2007; Shioya et al. 2008; Villar et al. 2008) 
except for points at z = 2.3 and 3.05 (Reddy et al. 2008) who used Ha luminosity function converted from 
UV luminosity function using a UV-Ha relation. 



